Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN M?eX style file vl.4) 



Acceleration time scale in an ultrarelativistic shock 



Janusz Bednarz 

Obserwatorium Astronomiczne, Uniwersytet Jagiellonski, ul. Orla 171, 30-244 Krakow, Poland 



1 February 2008 



o 
o 
o 

(N 



(N 
> 

O 

(N 

vn 
o 
o 
o 

Oh: 

6 



ABSTRACT 

The acceleration mechanism at uhrarelativistic shocks is investigated using the Monte 
Carlo simulations. We apply a method of discrete small amplitude particle momen- 
tum scattering to reproduce highly anisotropic conditions at the shock and carefully 
describe the acceleration mechanism. The obtained acceleration times equal 1.0 r^/c 
if the spectral index reach the value of 2.2, independent of physical conditions in the 
shock. Some other parameters of the acceleration process are also provided. 
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1 INTRODUCTION 

Observations carried out by the Burst and Transient Source 
Experiment show that GRBs originate from cosmological 
sources (Meegan et al. 1992 and Dermer 1992). Identification 
of the host galaxy for the GRB 971214 (Kulkarni et al. 1998) 
and several other bursts causes there is little doubt now 
that some, and most likely all GRBs are cosmological. These 
phenomena are surely related to ultrarelativistic shocks with 
the Lorentz factors 7 > 10^. 

Several papers suggested that ultrarelativistic shocks in 
GRBs could be sources of high energy cosmic rays (cf. Wax- 
man 1995, Vietri 1995), and simulations done by Bednarz & 
Ostrowski (1998) showed that such shocks are able to accel- 
erate charged particles and values of their energy spectral 
indices converge to a = 2.2 when 7 ^ 00 and/or magnetic 
turbulence amplitudes grow. Because the acceleration mech- 
anism is quite different from that in the non-relativistic and 
mildly relativistic regime we distinguish a class of ultrarela- 
tivistic shocks if their Lorentz factors 7 S> 1. 

Observations seem to confirm this mechanism. Waxman 

(1997) used a fireball model of GRBs and showed from the 
functional dependence of the flux on time and frequency that 
a = 2.3 ±0.1 in the afterglow of GRB 970228. Galama et al. 

(1998) made two independent measurements of the electron 
spectrum index in the afterglow of GRB 970508 which was 
very close to 2.2. 



2 ACCELERATION MECHANISM 

A particle crossing the shock to upstream medium has a mo- 
mentum vector nearly parallel to the shock normal. Then 
the particle momentum changes its inclination in two ways 
by: 1) scattering in an inhomogeneous magnetic field and 2) 
smooth variation in a homogeneous field component. Here- 



after, the mean defiection angle in these two cases will be 
denoted by AQs and Afi//, respectively. The first process is 
a diffusive one and the second depends on time linearly. That 
means that with increasing shock velocity, keeping other pa- 
rameters constant, Afls decreases slower as a square root of 
time in comparison with A^Ih- The Lorentz transformation 
shows that with 7^1 even a tiny angular deviation in 
the upstream plasma rest frame can lead to a large angu- 
lar deviation in the downstream plasma rest frame. Let us 
denote a particle phase by (p and the angle between mo- 
mentum and a magnetic field vector by 6 both measured 
in the downstream plasma rest frame. Values of these pa- 
rameters at the moment when a particle crosses the shock 
downstream determine if it is able to reach the shock again 
in the case of neglected magnetic field fluctuations down- 
stream of the shock. In fact a motion in the homogeneous 
magnetic field carries a particle in such a way that it can- 
not reach the shock again. The magnetic field fluctuations 
upstream of the shock perturbing the momentum direction 
lead to broadening the {cf>, 9) range that allows particles to 
reach the shock again. Thus, as we show below for efficient 
scattering, when AQ,h becomes unimportant in comparison 
to Afis, the spectral index and the acceleration time reach 
their asymptotic values. 

The discussed relation between AQ.H and AQ.s is repro- 
duced in our simulations and presented in Fig. 1. There are 
shown 11 points from 7 = 100 to 320 and three additional 
for 7 = 640, 1280, 2560. The expected linear dependence of 
these quantities can be noticed. 



3 NUMERICAL SIMULATIONS 

In simulations we follow the procedure used by Bednarz & 
Ostrowski (1996) with a hybrid approach used in Bednarz 
& Ostrowski (1998). Monoenergetic seed particles are in- 
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Figure 1. The relation between the mean deflection angle up- 
stream of the shock caused by the scattering in an inhomogeneous 
magnetic field (AHg) and by smooth variation in a homogeneous 
magnetic field (AO^f). Last three points for Afln below 1 • 10^^ 
represent 7 = 640, 1280, 2560 and yield a = 2.5, 2.3 and 2.2 
respectively. 
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shock (= particle density) is added to the time and mo- 
mentum bin depending on particle parameters, as measured 
in the shock normal rest frame. For the considered contin- 
uous injection after initial time, the energy cut-off of the 
formed spectrum shifts toward higher energies with time. 
The resulting spectra allows one to fit spectral indices and 
derive acceleration time in the shock normal rest frame in 
units of downstream rg/c ( rg - particle gyroradius in the 
homogeneous magnetic field component, c - speed of light; 
for details see Bednarz & Ostrowski 1996). We transform 
the acceleration time tacc to the downstream plasma rest 
frame. Hereafter, subscripts U or D mean that a parame- 
ter is measured in the upstream or downstream plasma rest 
frame respectively. We will use downstream rg as a distance 
and Tg/c as a time units. The magnetic field inclination to 
the shock normal upstream of the shock, ip, is measured in 
the upstream plasma rest frame. 

Let us denote the ratio of the cross-field diffusion co- 
efficient Kx to the parallel diffusion coefficient Ky as r (the 
value is measured in the plasma rest frame). Simulations 
prove that fluctuations upstream of the shock (measured by 
r^) and downstream of the shock (measured by r^) influ- 
ence the acceleration process independently. The minimum 
fluctuations upstream of the shock needed to run the accel- 
eration process efficiently tend to zero when 7 — > 00. We 
checked by simulations with different that its value does 
not influence the spectral index considerably for a given r^. 

Our scattering model is very simple but also universal. 
In the model we are not able to discuss gyroresonat scatter- 
ing. Upstream of the shock particles have not enough time to 
interact resonantly with low-frequency waves and the rough 
relation r - {5B/B f (cf. Blandford & Eichler 1987) cannot 
be deduce from the interaction there. However, for growing 
and fixed Ai9 the time between scattering acts decreases 
what is equivalent to increasing the magnetic field fluctua- 
tions. 



Figure 2. Simulated spectral indices as a function of magnetic 
field fluctuations upstream of the shock. Fluctuations downstream 
of the shock are neglected. The chosen value for 7 = 80 and 
ij) = 45° is pointed by a dashed line. A second-degree polynomial 
fit is also marked by a dashed line. 



jected at the shock and then their trajectories are derived 
in the perturbed magnetic field. The inhomogeneities are 
simulated by small amplitude particle momentum scatter- 
ing within a cone with angular opening Ai? less than the 
particle anisotropy ~ I/7 (cf. Ostrowski 1991). 

A particle is excluded from simulations if it escapes 
through the free-escape boundary placed far off the shock 
or reaches the energy larger than the assumed upper limit. 
These particles are replaced with the ones arising from split- 
ting the remaining high-weight particles, preserving their 
physical parameters. Particles that exist longer than the 
time upper limit for simulations are excluded from simu- 
lations without replacing. 

All computations are performed in the respective up- 
stream or downstream plasma rest frame. When particles 
cross the shock their parameters are transformed to the 
current plasma rest frame and the weighted contribution 
divided by the particle velocity component normal to the 



4 RESULTS 

In the following simulations we consider shocks with 7 = 
20, 40, 80, 160, 320, magnetic field inclinations V = 
15°, 30°, 45°, 60°, 75°, 90° and downstream values of mag- 
netic field fluctuations T° = 0, 1.0-10"^, l.l•10-^ 0.11, 0.69. 

Thus, as a flrst case we consider downstream conditions 
without magnetic field fluctuations. By simple data inspec- 
tion (cf. Fig. 2) we look for minimum where the spectral 
index reaches its limit of 2.2 and we apply this value in fur- 
ther simulations. The relation between . 7 and ip can be 
roughly fitted with the equation = 0.25 7^^'^ t/" in the 
considered range of shock parameters. We repeated sinm- 
lations for a number of cases with different 7 and V and 
r° ^ 0. The obtained results are in good agreement with 
the ones derived from the above equation up to — 0.11. 

Values of the acceleration time tacc for three amplitudes 
of magnetic field fiuctuations downstream of the shock are 
presented in Fig. 3. In the figure one can see the lack of 
change of tacc with ij), but it slowly decreases to the asymp- 
totic value with 7. In the simulations we have observed ten- 
dency of tacc to grow when a increases up to 2.3-2.4 and no 
further change if magnetic field fluctuations upstream of the 
shock grow. For < 0.11 the asymptotic value of the ac- 
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celeration time is close to rg/c. It occurs that rg/c is a good 
unit provided that the homogeneous magnetic field domi- 
nates the randomly component. Unfortunately, when this 
condition fails the meaning of tacc becomes unclear in the 
simulations then. For this reason we will not discuss further 
the case of = 0.69 any more. 

Approximate calculations of Gallant & Achterberg 
(1999) showed that t^/t§ ~ 1, where is the particle mean 
residence time upstream of the shock (upper index) as mea- 
sured in the upstream plasma rest frame (lower index), and 
D in stands for the downstream residence time. However, 
they were not able to consider the anisotropic particle mo- 
mentum distribution and our results in Fig. 4 transformed 
to the upstream plasma rest frame with t^/t^ within the 
range 0.01 — 0.1 are more adequate for real situations. Ad- 
ditionally, the above authors applied an extremely irregular 
magnetic field upstream of the shock represented by ran- 
domly oriented magnetic cells with field amplitude B and 
they measured time in the upstream unit of rg{B)/c. As a 
result they obtained that t^j/t^ could be much larger than 
1 in the case. 

Just before the spectral index reaches its minimal value 
(cf. Fig. 2) Af2s stabilizes near the limit which value does 
not further depend on the magnetic field inclination as is 
seen in Fig. 5. Momentum vectors of particles crossing down- 
stream of the shock have similar distributions as measured 
in the downstream plasma rest frame if A57s approaches 
the maximum value. Then, it follows that parameters we 
consider below depend only on r^. 



For growing r (r 



0, 1.0 • 10" 



1.1 • 10"^ 0.11) 



the acceleration time is constant and accompanied by a slow 
increase of the mean energy gain in one cycle downstream- 
upstream-downstream {AE/E)d = 0.89, 0.94, 1.0, 1.1, and 
a slight decrease of the fraction of particles that reach 
the shock again after crossing it downstream (An/n) — 
0.51, 0.50, 0.48, 0.44. Simultaneously the mean time a par- 
ticle spends downstream of the shock grows as = 
0.96, 1.0, 1.2, 1.35. Time that a particle spends upstream of 
the shock can be neglected in this rest frame as is visible in 
Fig. 4. It implies, approximately, tacc ~ tj^/{AE/E)D if one 
neglects correlations between these quantities (cf. Bednarz 
& Ostrowski 1996). Similarly we can roughly estimate the 
value of the energy spectral index of accelerated particles as 
cr~ l~ln{{An/n))/ln{{AE/E)D + l)- 

The simulated maximum distance in downstream 
medium the particle is able to depart from the shock and 
reach it again is, respectively, dj^ = 0.84, 1.5, 2.5, 4.0 (the 
values were derived from ~ 10^ events). We calculated 
the average values of sin^ 6 for returning particles wan- 
dering downstream of the shock and found, respectively, 
(sin^ e) = 0.676, 0.658, 0.650, 0.651. 
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Figure 3. The simulated acceleration time as a function of the 
shock Lorentz factor: a) without fluctuations downstream of the 
shock, b) with fluctuations downstream of the shock, c) fluctua- 
tions downstream of the shock dominate homogeneous magnetic 
field. Results for a given upstream magnetic field inclinations, 
given in a panel a), are joined with dashed lines. 
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5 DISCUSSION 

Because of newly found acceleration mechanism in ultrarel- 
ativistic shock waves we propose that some part of GRBs 
radiation could arise due to synchrotron radiation of elec- 
trons or electron pairs accelerated across the mechanism. 



Figure 4. The ratio of the mean time a particle spends upstream 
of the shock between shock crossings to the time it spends down- 
stream of the shock as a function of the shock Lorentz factor. 
It slightly decreases with growing upstream magnetic field in- 
clination. Dashed lines join points with a constant ■0. Apparent 
deviation of the point with 7 = 320, ip = 15° is real. 



* Below, we provide the respective series of simulated parameter 
for this sequence of 
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Figure 5. The mean deflection angle resulting from scattering 
in an inhomogeneous magnetic field upstream of the shock as a 

function of the shock Lorcntz factor 7 with chosen in a way as 
illustrated in Fig. 2. Six dashed lines for different ^ are identical. 
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We follow the internal shocks model of GRBs (cf. Kobayashi 
et al. 1997 for example). In the model two different shells 
have different Lorentz factors. The inner shell overtakes a 
slower outer shell and form a shock. The Lorentz factor of 
the shock as measured in the frame at rest with respect to 
the outer shell is assumed to be ~ 2. With only a part of ki- 
netic energy converted into the internal energy the particle 
energy distribution downstream of the shock will be non- 
thermal with, possibly, a substantial fraction of relativistic 
particles. As a result an amount of relativistic particles will 
be present in the shock. The particles can be accelerated 
across the mechanism presented in the paper and could be 
observed in the afterglow (cf. Waxman 1997, Galama et al. 
1998). 

We derived some parameters of the process that could 
be used in GRBs models. The acceleration time tacc = 
1.0 rg/c, measured in the downstream plasma rest frame in 
the unit of particle gyroradius in the homogeneous magnetic 
field component divided by the speed of light is the second 
important parameter besides the spectral index a — 2.2. The 
values of dj^ and t£ define the dimensions of the shock that 
allow the process to be effective, and the values of (sin^ 9) 
shows how the synchrotron radiation can influence the pro- 
cess. 

In the ejecta of the relativistic matter in the GRB model 
the outer shells can be faster than the following ones. In this 
case separated shocks with Lorentz factors reaching 7 ~ 10^ 
will be generated. The leading shock could produce seed pro- 
tons with energies of 10^* - 10^'' eV. These protons down- 
stream of the first shock can interact with the following one 
to be reflected with energy gains ~ 7^^ (cf. Gallant & Achter- 
berg 1999, Bednarz & Ostrowski 1999). For a constant re- 
flection probability the spectrum of these reflected highest 
energy particles, above 10^*^ eV, will be only the shifted in 
energy spectrum of seed particles with the universal spectral 
index ~ 2.2. 
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